高精度非线性格式WCNS的分析研究与其应用

首先将Fourier方法推广于高维方向研究了五阶精度WCNS的特性,并与其他高阶格式进行比较。分析结果表明WCNS的高精度特性普遍接近甚至好于迎风偏置五阶显式格式EUW-5与Pade′标准格式。然后开展了WCNS的应用研究,采用高效率的WCNS-E-5数值模拟了含强激波的高维复杂流场。算例包括二维高超声速边界层对自由流扰动波的吸收问题以及三维球头绕流问题。计算结果反映出WCNS-E-5对激波等间断的良好捕捉能力,图像清晰光滑,数据准确可靠。     

A numerical study of the motion of a wave running up a beach

The main difficulty for the numerical calculation of the wave running up a beach is the treatment of its moving water boundary. In this paper a scheme of turning the free boundary problem into a fixed boundary problem is designed. The calculated run-up height is consistent with the experiments. Some interesting wave phenomena are also found.     

类升力体外形高超声速粘性绕流拓扑结构分析

采用交替方向隐式分解的隐式NND格式求解全N-S方程模拟“类升力体”外形在高超声速下的大攻角流动,给出了“类升力体”外形表面极限流线随攻角变化的拓扑结构及40°攻角下垂直于体轴的横截面流线拓扑结构。结果表明:类升力体外形三维流场结构十分复杂,攻角从0°～40°变化时,背风面表面极限流线依次由不分离、开式分离向起始于鞍、结点组合的高阶奇点的分离方式转化,翼面横向分离亦随攻角增大而增大;垂直于体轴的横截面流动拓扑结构与文献[2]给出的理论分析一致。     

框架结构主动控制最优时滞研究

在结构振动主动控制中,时滞的大小对结构控制能量和响应均产生影响。在允许时滞的范围内,把结构控制能量和响应作为目标函数,数值仿真发现目标函数随着时滞的变化呈现曲线状态,这说明存在最优时滞。在目标函数处理方面,考虑到控制能量和响应数量级不一致,本文采用了无量纲化处理。同时,考虑到地震波的随机性,本文还推导了含控制力的改进的随机Newmark方法。最后给出了一个计算实例。实例表明本文所提出的方法是有效的。     

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

Mode-3 spontaneous crack propagation along functionally graded bimaterial interfaces

The effects of combining functionally graded materials (FGMs) of different inhomogeneous property gradients on the mode-3 propagation characteristics of an interfacial crack are numerically investigated. Spontaneous interfacial crack propagation simulations were performed using the newly developed spectral scheme. The numerical scheme derived and implemented in the present work can efficiently simulate planar crack propagation along functionally graded bimaterial interfaces. The material property inhomogeneity was assumed to be in the direction normal to the interface. Various bimaterial combinations were simulated by varying the material property inhomogeneity length scale. Our parametric study showed that the inclusion of a softening type FGM in the bimaterial system leads to a reduction in the fracture resistance indicated by the increase in crack propagation velocity and power absorbed. An opposite trend of increased fracture resistance was predicted when a hardening material was included in the bimaterial system. The cohesive tractions and crack opening displacements were altered due to the material property inhomogeneity, but the stresses ahead of the cohesive zone remained unaffected.     

A UNIFORMLY CONVERGENT DIFFERENCE SCHEME FOR THE SINGULAR PERTURBATION PROBLEM OF A HIGH ORDER ELLIPTIC DIFFERENTIAL EQUATION

ＡＵＮＩＦＯＲＭＬＹＣＯＮＶＥＲＧＥＮＴＤＩＦＦＥＲＥＮＣＥＳＣＨＥＭＥＦＯＲＴＨＥＳＩＮＧＵＬＡＲＰＥＲＴＵＲＢＡＴＩＯＮＰＲＯＢＬＥＭＯＦＡＨＩＧＨＯＲＤＥＲＥＬＬＩＰＴＩＣＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮ（刘国庆）（苏煜城）ＡＵＮＩＦＯ...     

Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD

Based on the successive iteration in the Taylor series expansion method,a three-point explicit compact difference scheme with arbitrary order of accuracy is derived in this paper.Numerical characteristics of the scheme are studied by the Fourier analysis. Unlike the conventional compact difference schemes which need to solve the equation to obtain the unknown derivatives in each node,the proposed scheme is explicit and can achieve arbitrary order of accuracy in space.Application examples for the convection- diffusion problem with a sharp front gradient and the typical lid-driven cavity flow are given.It is found that the proposed compact scheme is not only simple to implement and economical to use,but also is effective to simulate the convection-dominated problem and obtain high-order accurate solution in coarse grid systems.     

Fhree-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD

Based on the successive iteration in the Taylor series expansion method, a three-point explicit compact difference scheme with arbitrary order of accuracy is derived in this paper. Numerical characteristics of the scheme are studied by the Fourier analysisl Unlike the conventional compact difference schemes which need to solve the equation to obtain the unknown derivatives in each node, the proposed scheme is explicit and can achieve arbitrary order of accuracy in space. Application examples for the convectiondiffusion problem with a sharp front gradient and the typical lid-driven cavity flow are given. It is found that the proposed compact scheme is not only simple to implement and economical to use, but also is effective to simulate the convection-dominated problem and obtain high-order accurate solution in coarse grid systems.     

Some remarks about the resolution of high velocity flows near low densities

High velocity flows which are exposed to strong rarefaction waves and creating low densities regions in it present difficulties and inaccuracies for numerical resolution. In particular, variables related to the internal energy are wrongly evaluated. Use of classical schemes solving the Euler equations in conservative variables introduces significant errors in the determination of temperature. We recommend to employ a non-conservative formulation of the energy equation. Results found to be more accurate in using the present internal energy formulation. In order to have the formulation available for both shock and strong rarefaction waves, we propose a hybrid formulation of conservative and non-conservative ones, depending on a shock indicator. The results are compared with exact solutions and show a significant improvement of the accuracy. The method is then extended to two-dimensional cases. Received 28 March 1997 / Accepted 18 June 1997